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Abstract: We present high-precision numerical results for time-like Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi evolution in the MS factorisation scheme, for the first time up to 
next-to-next-to-leading order accuracy in quantum chromodynamics. First, we scrutinise 
the analytical expressions of the splitting functions available in the literature, in both x and 
N space, and check their mutual consistency. Second, we implement time-like evolution 
in two publicly available, entirely independent and conceptually different numerical codes, 
in x and N space respectively: the already existing APFEL code, which has been updated 
with time-like evolution, and the new MELA code, which has been specifically developed to 
perform the study in this work. Third, by means of a model for fragmentation functions, 
we provide results for the evolution in different factorisation schemes, for different ratios 
between renormalisation and factorisation scales and at different final scales. Our results 
are collected in the format of benchmark tables, which could be used as a reference for 
global determinations of fragmentation functions in the future. 
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1 Introduction 

In the framework of perturbative Quantum Chromodynamics (QCD), parton-to-hadron 
Fragmentation Functions (FFs) encode the information on how quarks and gluons are 
turned into hadrons [1, 2], Their precise knowledge is an essential ingredient in the quan¬ 
titative description of any hard-scattering process involving identified hadrons in the final 
state. Like Parton Distribution Functions (PDFs), FFs are non-perturbative quantities 
and, as such, they have to be determined from experimental data, typically in a global 
QCD analysis of a large variety of processes [3, 4], These analyses are all based on fac¬ 
torisation [5], which allows one to compute the relevant hard-scattering matrix elements 
perturbatively, and to absorbe the collinear singularities arising from the masslessness of 
partons into FFs. After factorisation, perturbative QCD corrections lead FFs to depend on 
the factorisation scale / ip , and this dependence obeys time-like Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) evolution equations [6-9]. 

Available experimental data to be included in a global QCD analysis of FFs span 
several orders of magnitude in energy. Beside rather old measurements (see ref. [4] for a 
review), new high-precision data are being produced copiously. On the one hand, they 
include multiplicities in fixed-target semi-inclusive deep-inelastic scattering (SIDIS) from 
HERMES [10] and COMPASS [11] experiments at 1.1 < Q 1 2 < 7.4 GeV 2 and 1.2 < Q 2 < 
22.4 GeV 2 respectively. On the other hand, they include production cross-sections of light 
hadrons in e + e _ collisions from BELLE [12] and BABAR [13] experiments at a center- 
of-mass energy y/s ~ 10.5 GeV, and in pp collisions from STAR [14] and PHENIX [15] 
experiments at y/s = 200 GeV and from CMS [16, 17] and ALICE [18] experiments up to 
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y/s = 7 TeV. Large Hadron Collider (LHC) data extend to unprecedented large values the 
energy reach, which will be even pushed forward by the future LHC Run-II. 

In order to determine FFs from these data sets, the evolution programs required for 
the multi-parameter global QCD analyses have to be numerically and conceptually under 
control. In principle, the demanded accuracy is the same as for PDFs: indeed, FFs and 
PDFs are on the same footing, and eventually they could be determined simultaneously in 
a fit to experimental data. This may be of particular interest in the case of spin-dependent 
PDFs, since a large amount of the experimental information used for their determination 
comes from SIDIS data with longitudinally polarised beams and targets. In this case, 
the potential cross-talk between spin-dependent PDFs and FFs could then be profitably 
addressed in a fit of both these quantities based on a mutually consistent methodology. 

So far, high-precision numerical results for time-like evolution have not been presented 
in a systematic way. In very much the same spirit of previous studies for space-like evolu¬ 
tion [19, 20], the goal of this paper is to fill this gap: specifically, we present results obtained 
in the MS factorisation scheme and, for the first time, up to next-to-next-to-leading order 
(NNLO) accuracy in QCD. Our study aims at providing a reference for future global QCD 
analyses of FFs. 

Our goal is achieved in three steps. First, we check the mutual correspondence of 
time-like splitting functions in the literature in both x and N space whenever available. 
Second, we implement them in two entirely independent and conceptually different evolu¬ 
tion programs: APFEL (A PDF Evolution Library) [21, 22], an already existing evolution 
package which in version 2.3.0 has been updated with time-like evolution, and MELA (Mellin 
Evolution LibrAry), a new program which has been developed specifically for the purpose 
of this paper. Third, we use a model for FFs with a sufficient degree of realistic behavior in 
order to obtain our reference results. Provided perfectly controlled conditions, obtaining 
the same results irrespective of the procedure followed in the two codes provides a strong 
check of the correctness of both our implementations and our results. 

Note that time-like evolution is performed in the MS factorisation scheme by two 
publicly available programs, QCDnum [23] and ffevol [24], and by private codes used for re¬ 
cent global determinations of FFs [25-27], though only up to next-to-leading order (NLO) 
accuracy. Furthermore, these codes differ among each other for physical and technical 
assumptions, hence a comparison between the results obtained either with QCDnum and 
ffevol or with one of these programs and published parameterisations is not straight¬ 
forward (or even not possible). In this respect, our results provide a reference for future 
comparisons with these codes, calling for a dedicated effort beyond the scope of this work. 

The paper is organised as follows. In section 2, we review the structure of the DGLAP 
evolution equations, and we scrutinise the expressions of the time-like splitting functions 
available in the literature. In section 3, we describe the numerical implementation of the 
time-like evolution in MELA, we discuss the setup conditions for our benchmark versus APFEL, 
and we present corresponding results and accuracy. Two appendices complete our paper. 
In appendix A, we collect the numerical results of our study in the format of benchmark 
tables. Finally, in appendix B, we provide a minimal set of instructions for downloading 
and running the MELA benchmark code. 
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2 Time-like evolution 


The evolution of parton-to-hadron FFs is described by 2n/ + 1 DGLAP equations 

3^2 - d H x ^f) = T, f x ~ P f (2b a s(MF)) D j > ( 2i1 ) 

where rif is the number of active flavours, the index j runs over partons, D ,■* is the FF 
for the parton i to fragment into a hadron h, x is the scaled energy of the hadron h (i.e. 
the fraction of the parton four-momentum taken by the hadron h), y f is the factorisation 
scale, a s is the QCD running coupling, and Pji are the time-like splitting functions. These 
allow for a perturbative expansion of the form: 

P J1 {y,a s ) = Y J O k s +l P ( £\y), ( 2 . 2 ) 

fc =0 

where we have defined a s = a s /(47r). From considerations based on charge conjugation 
and flavour symmetry, eqs. (2.1) are usually rewritten into 2re/ — 1 equations 

<91n/r 2 P ns-,±,v( x i Pf) = T ±,v (x, y 2 F ) <g> D^ s . ± v (x, y F ) (2.3) 

describing the independent evolution of non-singlet quark FF combinations, = {D^z t 

DP) — (Dg ± Dg-.) and s . v = 'Y^L 1 {Dq i — D f r j .), and a system of two coupled equations 

_d_(D*(x,&) W P™ 2 n f P*n\ (D*(x,tf,)\ 

dlnyl [d^x^D ) {&jP « P« ) ® \D" g {x^l) ) { ■ ] 

describing the evolution of gluon and singlet fragmentation functions, D g and D ^ = 
Y2*i=i {Pq, + P<jJ- The shorthand notation <8) stands for the convolution product 

/ (x) <S> g(x) = [ —f(y)g(~'] ■ (2.5) 

Jx y \yj 

The solution of eqs. (2.1) can be performed either in x or in N space; in N space, 
the system of integro-differential equations (2.1) becomes a system of Ordinary Differential 
Equations (ODEs) of the form 

^-^D^{N,y 2 F ) = Y J Pqi (IV,«,(/£)) P) {N,&) . (2.6) 


The splitting functions in the two spaces can be related to each other via the Mellin 
transform 


and its inverse 


Pji{N ,a s )= / dyy 


/' 

Jo 


..TV -1 


Pjiilli a s) 


1 

2ni 


rc-t 
J c—v 


c+ioo 


Pji(y,a s ), Ne C, 

(2.7) 

dN y~ N Pji(N, a s ), 

(2.8) 
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where the real number c has to lie to the right of the rightmost singularity of P Jt in the 
complex plane. Fragmentation functions in x and A space are related with each other by 
the same transformations as in eqs. (2.7)-(2.8). 

At leading order (LO), time-like splitting functions are identical to their space-like 
counterparts, while they differ at higher orders. At NLO, explicit expressions for the 
complete set of time-like splitting functions in eqs. (2.3)-(2.4) are collected in refs. [28, 29] 
(x space), in ref. [30] (A space) and ref. [31] (both x and A space), though with rather 
different notations. 1 All results are given in the MS factorisation scheme. 

Time-like splitting functions have been computed up to O(o.g), i.e. NNLO accuracy, 
always in the MS factorisation scheme in refs. [35-37]. An uncertainty still remains on the 
exact form of P r \ g . Indeed, this was determined by means of a relation between known 
NLO evolution kernels for photon- and Higgs-exchange structure functions in deep-inelastic 
scattering, and their counterparts in semi-inclusive annihilation [32, 38], supplemented with 
constraints arising from momentum sum rule and supersymmetric relations for the choice 
Ca = Cf = nf of colour factors. The latter fix the form of Pq g J except for the offset 
quantified by Eq. (38) in ref. [37], which does not affect the logarithmic behaviour of Pq g ' 
at small and large momentum fractions [37] and consequently the validity of our study. Note 
finally that coefficient functions are known at NNLO only for e+e~ annihilation [31, 39-41]: 
this will thus limit the potential of a global determination of fragmentation functions at 
NNLO. 

For doubt’s sake, we checked that all the aforementioned results, in both x and A space, 
fully agree with each other up to NNLO accuracy. At LO, such a check is straightforward 
due to the extreme simplicity of the expressions involved, and we found perfect agreement 
among all the results considered. At NLO and NNLO, instead, time-like splitting functions 
are much more complicated than at LO. Also, at NLO the check is complicated by the 
fact that different notations, not directly comparable, and different FF basis are used in 
refs. [29-31]. Be that as it may, we considered the basis of refs. [31, 35-37]: at NLO, this 
is {P+,P~,PP S ,P^,P^,P^}^ at NNLO, this is {P+, P~, P v , P ps , P q s, P§s}. 2 This 
corresponds to the basis entering eqs. (2.3)-(2.4), provided that P qq = P+ + P ps . We refer 
to chapter 4 of ref. [29] for another example of FF basis. 

At NLO, we carried out the comparison in two steps. First, we checked that the x-space 
results in refs. [29, 31] and the A-space results in refs. [30, 31] are identical. We performed 
this check both analytically and numerically. On the one hand, in order to deal with the 
analytic notation in ref. [31], we have used the definitions of harmonic polylogarithms and 
harmonic sums provided in refs. [42, 43]. On the other hand, we have used the package 
presented in ref. [44] to handle harmonic sums numerically. Second, we checked that the 
x-space results in ref. [29] correspond to the A-space results in ref. [30], and that x- and 
A-space results in ref. [31] correspond to each other. We performed this check numerically, 

1 Results in ref. [29] were obtained for the first time in ref. [28]. Well-known misprints in ref. [28] have 
been amended in ref. [29], see also refs. [32, 33]. Note that A-space results in ref. [30] were obtained by 
applying the Mellin transform, eq. (2.7), to the x-space results in ref. [28], while those in ref. [31] were 
derived directly in A space using the method presented in ref. [34], 

2 Note that, at NLO, there are only six independent splitting functions because P v = P~. 
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Mitov-Moch x-space vs. Ellis-Stirling-Webber (NLO) 


Mitov-Moch vs. Glueck-Reja-Vogt (NLO) 
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Figure 1: (Left panel). The relative difference A P/P between the x-space splitting functions of 
ref. [31] and those of ref. [29]. (Right panel). The absolute relative difference |AP/P| between the 
IV-space splitting functions of ref. [31] and those of ref. [30]. 


by transforming the IV-space expressions in x space, and by comparing the results with 
the corresponding x-space expressions (for details about the implementation of the inverse 
Mellin transform, see section 3.1 below). 

At NNLO, we checked that the x- and N -space expressions provided in refs. [35-37], 
and available in ref. [45] in the format of Fortran subroutines, correspond to each other. 
Again, we performed this check numerically, by transforming the IV-space expressions in 
x space, and by comparing the results with the corresponding x-space expressions. Note 
that we considered the approximate representation (or parameterisation) of NNLO time¬ 
like splitting functions provided in ref. [45], consistently in x and N space. Indeed, the 
exact expressions of the NNLO splitting functions are rather complex and these will lead to 
very lengthy computations once implemented in a code for numerical evolution like APFEL. 
It was checked in ref. [37] that, except for values of x very close to zeros of the splitting 
functions, such approximate expressions deviate from the exact ones by less than one part 
in a thousand. 

In the left panel of figure 1, we show the relative difference between the x-space splitting 
functions of ref. [29] and those of ref. [31] over a sensible range in x. In the right panel of 
figure 1, we show the absolute value of the relative difference between the IV-space splitting 
functions of ref. [30] and those of ref. [31] over a wide portion of the complex plane. In the 
left panel of figure 2, we show the relative difference between the exact x-space splitting 
functions in ref. [29] and the numerical inversion of the IV-space expressions of ref. [30] (at 
NLO). In the central panel of figure 2, we show the same comparison but for the splitting 
functions in ref. [31] (at NLO). In the right panel of figure 2, we show the same comparison, 
but for the NNLO splitting functions in refs. [35-37, 45]. 

The results displayed in figures 1 and 2 allow us to draw the following conclusions. 

• The agreement between the expressions of time-like splitting functions at NLO, in x 
and N space separately, is optimal. In the x-space case (left plot in figure 1), the 
relative difference between them is A P/P ~ 10 -10 . In the IV-space case (right plot 
in figure 1), the range for the absolute relative difference between them is 10 -11 < 
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Glueck-Reya-Vogt vs. Ellis-Stirling-Webber (NLO) 


Mitov-Moch, x-space vs. N-space (NLO) 


Almasy-Moch-Vogt, x-space vs. N-space (NNLO) 



Figure 2: (Left panel). The relative difference A P/P between the x-space splitting functions of 
ref. [29] and the inverse Mellin transform of the TV-space splitting functions from ref. [30]. (Central 
panel). The relative difference A P/P between the x-space and TV-space splitting functions of 
ref. [31] (at NLO). (Right panel) The relative difference A P/P between the x-space and iV-space 
splitting functions of refs. [35-37, 45] (at NNLO). 


\AP/P\ < 10 -5 . The fact that the values of \AP/P\ in the JV-space case cover a 
larger spread than the values of A P/P in the x-space case is a consequence of the 
numerical evaluation of harmonic sums with the package provided in ref. [44]. Indeed, 
this is a further source of numerical uncertainty, which, however, is well under control. 

• The agreement between the inverse Mellin transform of N-space time-like splitting 
functions and their x-space counterparts is also good, both at NLO and NNLO. 
Indeed, in figure 2 the relative difference between them is A P/P < ICO 4 , irrespective 
of the splitting function and the perturbative order. Note however that the value of 
A P/P is larger in figure 2 than in figure 1. In the former case, there is an additional 
uncertainty related to the numerical inverse Mellin transform of IV-space splitting 
functions back to x space. 

We conclude that the expressions of the splitting functions available in the literature 
are perfectly consistent with each other. 3 

3 Numerical implementation and benchmark 

In this section, we present our implementation of time-like evolution in the MS factorisation 
scheme up to NNLO accuracy. The discussion is organised in three steps. First, we 
briefly discuss the numerical solution of DGLAP equations and their implementation in 
two different programs. Second, we define the setup conditions. Third, we present the 
results of a benchmark between our two codes. 

3 We would like to draw the reader’s attention on a couple of minor misprints that we came across in the 
literature during the checks we performed. In the arXiv version of ref. [42], there should be a minus sign, 
rather than a plus sign, in front of Lbfx) in the definition of H( 1,0; x), eq. (11). In ref. [44], in the definition 
of S-k(N), second relation in eq. (10), the squared bracket should be closed before £(fc), rather than after, 
see also eq. (46) of ref. [46]. Finally, a couple of misprints affecting both x- and N-space expressions of P qs 
in ref. [31], eqs. (C.10) and (B.10) respectively, have been corrected in a revised version recently submitted 
to the arXiv. 
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3.1 The codes 


We have implemented the time-like evolution in two entirely independent and conceptually 
different programs, based respectively on the solution of DGLAP equations (2.1)-(2.6) in 
x and N space. Provided perfectly controlled conditions, our goal is to obtain the same 
results irrespective of the procedure followed in the two codes. This provides a strong check 
of the correctness of both our implementations and our results. 

As far as the x-space solution is concerned, we have implemented it in the already 
existing APFEL library [21]. Specifically, APFEL provides a framework in which the DGLAP 
equations are solved in x space by means of higher-order interpolation techniques, followed 
by the Runge-Kutta solution of the resulting discretised evolution equations. We refer the 
reader to ref. [21] for technical details on the implementation. 

As far as the N -space solution is concerned, we have developed a new dedicated pro¬ 
gram: MELA, Mellin Evolution LibrAry. MELA exploits the fact that, in N space, the integro- 
differential DGLAP equations (2.1) are turned into the more simple ODEs (2.6). While 
eqs. (2.1) can be solved only by means of numerical methods, eqs. (2.6) allow for a sim¬ 
ple analytical solution [47]. The numerical evaluation of the latter is immediate and, in 
principle, infinitely accurate. Of course, one should perform the numerical inversion of 
the solution from N space back to x space, but this is technically much easier than solv¬ 
ing eq. (2.1) directly. In MELA, the inverse Mellin transform is performed by means of a 
numerical implementation of eq. (2.8) based on the Talbot-path algorithm [48]. 

Note that the solution of eqs. (2.6) requires the knowledge of the analytical expressions 
of the Mellin moments of the FFs over the whole complex plane. Hence, the parameter- 
isation of FFs, usually given in x space, should be sufficiently simple to allow for an 
analytical Mellin trasform. Of course, this greatly restricts the range of application of the 
N-space solution of the DGLAP equations. In order to bypass this limitation, the so-called 
FastKernel method [49], which allows one to evolve x-space distributions using the N- 
space approach by means of interpolation techniques, has been developed. This method 
has been implemented in MELA for completeness, even though this is not required by our 
study. Further investigations based on the extension of the FastKernel method to the 
time-like case will be left for future studies. 

3.2 The setup 

Our goal is to consistently compare APFEL and MELA. In order to do that, we need to choose 
a common setup for the evolution. Specifically, we use the following value of the strong 
coupling at the charm mass m c : 

a s (m c ) = 0.35, (3.1) 

and the following values for the heavy quark masses: 

m c = 1.43 GeV, mj, = 4.3 GeV. (3-2) 

The top quark mass mt does not need to be specified because the top-mass threshold will 
never be crossed in our computations. 
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We take the parameterisation for the initial-scale FFs from ref. [25]. In particular, we 
consider the 7r + FFs at NLO, eqs. (14)-(16), with the parameter values collected in table VI 
of ref. [25] 

D?(x,f4) = - x) 1 - 370 , 


D l + ( x ,Vo) = D^ + (x,n g) = Df(x,nl) 

Df(x,nl) = AjV' 943 (l-x) 8 . 


Df{x^l)=N t 


2\ _ A7'-7r+ „0.718 


X 


*(1 — x ) 


6.266 


(3.3) 


Here /Xq = 1 GeV 2 and the normalisation constants N£ + , VJ + and Ng + are such that 

f dx xD^ + {x, = 0.401, [ dx xD^ + (x, = 0.094, f dx xD^ + (x, = 0.238. 

Jo Jo Jo 

(3.4) 

Note that we consider only gluon and light quark FFs: indeed, in our programs heavy- 
quark components of FFs are dynamically generated at the corresponding thresholds. For 
this reason, they do not need to be parameterised. Consistently, we include matching 
conditions in the treatment of flavour threshold crossing in the evolution (2.1): in this 
respect, we note that, differently from PDFs, FFs undergo a non-zero matching at the 
heavy-quark thresholds already at NLO [50], while the matching at NNLO is presently 
unknown. 

We would like to stress that the choice of the parameterisation given by eq. (3.3) is 
arbitrary. Essentially, this is motivated by its simplicity, which allowed us to easily obtain 
the corresponding analytic Mellin trasform required by MELA. Therefore, this parameterisa¬ 
tion should not be considered more reliable that any other, and the FFs given in eq. (3.3) 
should be regarded just as a set of test functions with some degree of realistic behaviour. 

3.3 The results 

The benchmark between APFEL and MELA is performed by comparing the result of the 
evolution of the set of FFs (3.3) for two different factorisation schemes, for different ratios 
between renormalisation and factorisation scales and at different final scales. 

As far as the factorisation scheme is concerned, at NLO we consider two options: 
the Fixed-Flavour-Number Scheme (FFNS) with rif = 3, in which the number of active 
flavours remains fixed and no heavy quark FF is generated during the evolution, and the 
Variable-Flavour-Number Scheme (VFNS), in which the heavy-quark FFs are dynamically 
generated as the evolution scale crosses the corresponding heavy-quark thresholds. The 
lack of knowledge on matching conditions prevents us to provide results in the VFNS at 
NNLO: hence, in this case, we perform the evolution only in the FFNS. 

As a further cross-check, we also consider the more general case in which the factori¬ 
sation scale at which FFs are evaluated, and the renormalisation scale /xr, at which a s 
is evaluated, take different values {fip / Mr)- In this case, the effective splitting functions 
depend on the coefficients of the perturbative expansion of the QCD /3-function and on 
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APFEL vs. MELA: FFNS at NLO 
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Figure 3: The ratio -Rmela/Apfei/z) = -^iMELA^V-^APFEL^) f° r the comparison between 
the positive charged pion fragmentation functions of gluon, up, strange and charm quarks evolved 
at different values of y? F with APFEL and MELA. Results refer to different values of the ratio of the 
renormalisation and factorisation scales, in the FFNS at NLO. 


powers of In (isp/n\) [47]. In addition, if /zp / hr, the matching of a s is no longer per¬ 
formed at the heavy-quark threshold, thus giving rise to discontinuities in the evolution of 
a s . In our benchmark we consider three different values for the ratio Hr/ l^p'- 1/2,1, 2. 

We evolve the set of FFs, eq. (3.3), from /// = 1 GeV 2 to four different final scales, 
namely Hf = 10, 100, 1000, 10000 GeV 2 . As for the momentum fraction, we look at 
the range x > 0.01, consistently with the kinematic cut usually imposed in global QCD 
analyses of FFs. Actually, beyond LO, time-like splitting functions have a much more 
singular behaviour than their space-like counterparts as x —> 0. Specifically, they present a 
double-logarithm enhancement with leading terms of the form a™ ln 2n_2 x (corresponding 
to poles of the form a™(N — l) 1 ” 2 ™ in N space). Despite large cancellations between leading 
and non-leading logarithms at non-asymptotic values of x, the resulting small-x rise in the 
time-like case dwarfs that of the space-like case. This justifies the rather restricted range 
of momentum fractions we look at. 

In figure 3, we show the ratio between gluon, up quark and strange quark FFs evolved 
at NLO with APFEL and MELA from Ho = 1 GeV 2 to Hf = 10, 100, 1000, 10000 GeV 2 
for three different values of the ratio n\/ ii 2 f in the FFNS. In figure 4, we show the same 
comparison as in figure 3, but in the VFNS; here we display also the ratio for the charm 
quark FF, which is dynamically generated at the charm threshold. In figure 5, we show 
the same quantity as in figure 3, but at NNLO. All results are in the MS scheme. 

Figures 3-4-5 show that the agreement between APFEL and MELA is optimal over all 
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APFEL vs. MELA: VFNS at NLO 
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Figure 4: Same as figure 3, but in the VFNS. 
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Figure 5: Same as figure 3, but at NNLO. 


the (x, ll 2 f ) region explored, the accuracy being below the 10 -4 level. A slight worsening 
is observed only for the strange and charm quark FFs in the very large-x region. Here the 
absolute values of these distributions become extremely small and numerical effects start 







































coming in. 

We conclude that proper implementations of both x- and iV-space methods for the 
solution of the time-like DGLAP equations have been achieved in APFEL and MELA up to 
NNLO. The reliability of our results has been carefully checked step by step in order to 
exclude any inconsistency from the core of the evolution (splitting functions) up to the 
complete evolution mechanism. In order to facilitate a possible comparison with other 
codes, in appendix A we present a collection of tables reporting the numerical values of 
FFs at some reference values of x and /ip. The tables in appendix A, as well as the plots 
in figures 3-4-5, can be reproduced by running the benchmark suite of MELA, as explained 
in appendix B. 

4 Conclusions and outlook 

We have presented numerical results of the time-like evolution of FFs in the MS factorisa¬ 
tion scheme for the first time up to NNLO accuracy in QCD. A high-precision comparison 
between x- and N-space solutions of the DGLAP evolution equations has been provided 
based on two entirely independent and conceptually different programs: APFEL, which has 
been updated for handling tinre-like evolution, and MELA, which has been specifically de¬ 
veloped for the purpose of this paper. 

This study has made us scrutinise the expressions of the time-like splitting functions 
available in the literature, in both x and N space. We have checked that all available 
results in the literature are mutually consistent. 

We have obtained an excellent numerical agreement between APFEL and MELA results, 
achieving a relative accuracy well below 10 -4 . The stability and reliability of the evolution 
codes has also been extensively tested upon various theoretical assumptions, including 
fixed- and variable-flavour number schemes and different factorisation and normalisation 
scale ratios. Above all, we managed to control the evolution of FFs at a level of accuracy 
which is comparable to that demanded for PDFs. 

Our results aims at providing a reference for future global analyses of FFs; these may 
include a determination based on the NNPDF methodology [51, 52], thanks to the efficiency 
and flexibility of APFEL. In the future, we hope our results will be compared with other 
available codes, namely QCDnum and ffevol: this will call for a dedicated effort, which 
should benefit from the collaboration of the authors of these programs. 


APFEL is available from its HepForge website: 

http://apfel.hepforge.org/ 

MELA is available from: 


http://apfel.hepforge.org/mela.html 
We provide basic instructions for downloading and running it in appendix B. 
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A Numerical Results 

In this appendix, we present a collection of tables with the values of the FFs evolved as 
explained in section 3. We recall that we use the following values for the strong coupling 
at the charm mass m c and the charm and bottom quark masses: 

<Vs(m c ) = 0.35 , m c = 1.43 GeV , m& = 4.3GeV. (A.l) 

The parameterisatons for the initial-scale FFs are taken from ref. [25] with the parameter 
values collected in table VI in that reference. 

The following tables are meant to be used for future comparisons with other computa¬ 
tions. The values contained in the tables, as well as the comparison plots in figures 3-4-5, 
should be reproducible by running the MELA benchmark code as explained in appendix B. 

B The MELA evolution code 

MELA (Mellin Evolution LibrAry) is a A-space evolution program developed specifically 
for the computation presented in this paper. The aim of MELA is to provide a simple and 
user-friendly cross-check of the time-like x-space evolution available in APFEL since v2.3.0. 

Download 

The last stable release is MELA 1.0.0 and it is available from the HepForge website: 

http://apfel.hepforge.org/mela.html 

In order to download the latest version and decompress the code locally, open a terminal 
and execute the following commands: 

wget http://apfel.hepforge.org/downloads/MELA-1.0.0.tar.gz 
tar xvf MELA-1.0.0.tar.gz 


Running the benchmark code 

In order to produce the plots presented in section 3.3, as well as the benchmark tables 
presented in appendix A, execute the following commands: 

cd ./MELA/usr/ 

./Benchmark.sh 
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X 

Df{x 

,Qo ) 

Df{x 

,Qo ) 

Df(x 

,Qo ) 

1.0- 

l 0 -°2 

1.4028 • 

10- 01 

7.5485 • 

10-oi 

-6.7003 

- IO -02 

5.0- 

IQ-02 

1.0215 • 

lO+oo 

1.2618 • 

lO+oo 

5.6882 - 

10-oi 

1.0- 

io- 01 

9.6274 • 

10- 01 

1.0816 ■ 

10+oo 

5.4671 ■ 

10-oi 

2.0- 

10-oi 

7.2835 • 

IO” 01 

7.4260 - 

10-oi 

3.8267 • 

10-oi 

3.0- 

io- 01 

4.0737 • 

IO- 01 

5.1131 - 

10-oi 

2.1467- 

10-oi 

4.0- 

10-oi 

1.6711 • 

10~ 01 

3.5699 ■ 

10-oi 

9.8314- 

Iq-02 

5.0- 

io- 01 

4.9557 • 

l 0 -°2 

2.4692 • 

10-oi 

3.5922 - 

Iq-02 

6.0- 

IO" 01 

1.0252 • 

IQ-02 

1.6227- 

10-oi 

9.7576 • 

l 0 -° 3 

7.0- 

io- 01 

1.5593- 

CO 

O 

1 

O 

9.5732 • 

Iq-02 

1.6906 - 

l 0 -° 3 

8.0- 

IO" 01 

2.6359 • 

o 

1 

O 

4.5769 • 

Iq-02 

1.3035 • 

l 0 -° 4 

9.0- 

10-oi 

3.9247 • 

10-° 5 

1.3028 • 

Iq-02 

1.3262 • 

Iq-06 


Table 1: Numerical values of gluon, up and strange FFs evolved at NLO in the FFNS up 
to /4 = 10 GeV 2 , = 1/2, a s (n 2 F ) = 0.2770096. 


X 

Dg + (x,Qo) 

Df{x 

,Qo ) 

Df{x 

, Q o) 

1.0- 

Iq-02 

4.9764 • 

10-oi 

9.7894 • 

10-oi 

1.5823 • 

10-oi 

5.0- 

Iq-02 

8.6824 • 

IO- 01 

1.1681 • 

lO+oo 

4.7168 • 

10-oi 

1.0- 

10-oi 

8.3316- 

IO- 01 

1.0034 • 

10+oo 

4.6351 ■ 

10-oi 

2.0- 

IO- 01 

6.7795 • 

10-oi 

7.1424 • 

10-oi 

3.5073 • 

10-oi 

3.0- 

10-oi 

4.0174 • 

10-oi 

5.0877 • 

10-oi 

2.0860 - 

10-oi 

4.0- 

IQ-Oi 

1.7301 ■ 

IQ-Oi 

3.6378 • 

10-oi 

9.9818- 

IQ-02 

5.0- 

10- Ql 

5.3680 - 

Iq-02 

2.5559 • 

10-oi 

3.7734 • 

Iq-02 

6.0- 

IQ- 01 

1.1609 ■ 

iq-02 

1.7010- 

10-oi 

1.0539 - 

l 0 -° 2 

7.0- 

IQ- 01 

1.8188- 

1Q-0S 

1.0175 • 

10-oi 

1.8733 • 

l 0 -° 3 

8.0- 

IQ-Oi 

2.9568 • 

10 04 

4.9561 - 

l 0 -° 2 

1.4898 • 

l 0 -° 4 

9.0- 

IQ- 01 

4.1078- 

IQ- 05 

1.4556 • 

l 0 -° 2 

1.6551 - 

Iq-06 


Table 2: Same as table 1: FFNS, fi 2 F = 10 GeV 2 , [i 2 R jfj? F = 1, a s {^ 2 F ) = 0.2393111. 


X 

Dg + (x,Qo) 

Df{x 

! Qo ) 

Df{x 

,Qo ) 

1.0- 

Iq-02 

5.8087 - 

10-oi 

1.0412 • 

lO+oo 

2.2016 • 

10-oi 

5.0- 

Iq-02 

7.7199 - 

10-oi 

1.1158 • 

10 +oo 

4.1428 • 

10-oi 

1.0- 

10- Ql 

7.6289 ■ 

10-oi 

9.6544 • 

IO- 01 

4.1913 • 

10-oi 

2.0- 

10- Ql 

6.6472 • 

10-oi 

7.0473 • 

IO- 01 

3.3668 • 

10-oi 

3.0- 

10- Ql 

4.1434 • 

10-oi 

5.1305 • 

IO- 01 

2.0858 • 

10-oi 

4.0- 

IQ-Oi 

1.8581 ■ 

IO" 01 

3.7269 • 

io- 01 

1.0285 • 

10-oi 

5.0- 

10-oi 

5.9703 - 

Iq-02 

2.6492 • 

IO- 01 

3.9815 • 

Iq-02 

6.0- 

10- 01 

1.3277 - 

Iq-02 

1.7818 • 

IO- 01 

1.1349- 

iq-02 

7.0- 

IQ-Oi 

2.0841 ■ 

1 0 03 

1.0788 • 

IO- 01 

2.0572 • 

IQ- 03 

8.0- 

IQ- 01 

3.2211 ■ 

IO- 0 4 

5.3420 ■ 

IQ-02 

1.6750 - 

l 0 -° 4 

9.0- 

IQ- 01 

4.3275 • 

10-os 

1.6124 • 

I— 1 

o 

1 

o 

to 

1.9614 • 

IQ-06 


Table 3: Same as table 1: FFNS, [i 2 F = 10 GeV 2 , [i 2 R /[j? F = 2, a s {^ 2 F ) = 0.2110517. 
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X 

Dg + (x,Q 0 ) 

Df (x 

3 Go) 

Df(x 

,Go) 

1.0- 

10-°2 

1.4340 ■ 

10 +oo 

1.3327 • 

lO+oo 

5.2393 • 

10-oi 

5.0- 

Iq-02 

1.4491 ■ 

lO+oo 

1.4809 ■ 

10 +oo 

8.3349 ■ 

10-oi 

1.0- 

nr 01 

1.0512 ■ 

lO+oo 

1.1286 ■ 

10 +oo 

6.4251 ■ 

10-oi 

2.0- 

to- 01 

5.5687 • 

10-oi 

6.8324 • 

IO- 01 

3.5879 - 

10-oi 

3.0- 

to- 01 

2.4010 ■ 

io- 01 

4.3584 • 

IO- 01 

1.7432 • 

10-oi 

4.0- 

10-oi 

8.0347 • 

l 0 -° 2 

2.8747 • 

IO- 01 

7.1450- 

Iq-02 

5.0- 

to- 01 

2.0548 • 

l 0 -° 2 

1.8837 • 

IO- 01 

2.3713- 

Iq-02 

6.0- 

10-oi 

4.1359- 

1Q-0S 

1.1647- 

IO- 01 

5.8741 - 

l 0 -° 3 

7.0- 

to- 01 

8.0031 - 

10“ 04 

6.3757 • 

IQ- 02 

9.2263 • 

10” 04 

8.0- 

to- 01 

1.7251 ■ 

IO" 04 

2.7544 • 

IQ- 02 

6.3248 • 

IQ-OS 

9.0- 

to- 01 

2.0689 ■ 

10 -os 

6.6359 • 

CO 

o 

1 

O 

5.4430 - 

10-07 


Table 4: Same as table 1: FFNS, / i 2 F = 10 2 GeV 2 , n 2 R /n 2 F = 1/2, a s ([i 2 F ) = 0.1829054. 


X 

Df(x,Qo) 

Df{x 

3 Go) 

Df{x 

, Go) 

1.0- 

Iq-02 

1.5369 • 

10+ 00 

1.4617- 

lO+oo 

6.5377 - 

10-oi 

5.0- 

Iq-02 

1.2399 • 

lO+oo 

1.3476 ■ 

10 +oo 

6.9561 ■ 

IO - 01 

1.0- 

IO - 01 

9.2734 • 

10-oi 

1.0438 • 

10 +oo 

5.5200 ■ 

IO - 01 

2.0- 

10-oi 

5.2786 • 

10-oi 

6.6130- 

IO - 01 

3.3263 • 

IO - 01 

3.0- 

IO - 01 

2.4154 • 

10-oi 

4.3762 • 

IO - 01 

1.7148 • 

IO - 01 

4.0- 

IQ-Oi 

8.5126- 

Iq-02 

2.9581 - 

10-oi 

7.3565 • 

Iq-02 

5.0- 

IO - 01 

2.2822 • 

Iq-02 

1.9707 • 

IO - 01 

2.5309 • 

Iq-02 

6.0- 

IO” 01 

4.7668 • 

Iq-OS 

1.2362 • 

IO - 01 

6.4603 - 

Iq-OS 

7.0- 

IQ - 01 

9.2483 • 

10 04 

6.8785 • 

l 0 -° 2 

1.0440 ■ 

IQ - 03 

8.0- 

io- 01 

1.9406 • 

10 04 

3.0385 • 

l 0 -° 2 

7.4134 • 

10 -os 

9.0- 

IQ - 01 

2.2911 ■ 

10 -os 

7.5963 • 

CO 

o 

1 

O 

7.0588 • 

Iq-07 


Table 5: Same as table 1: FFNS, [i 2 F = 10 2 GeV 2 , fj? R /n F = 1, a s {/j? F ) = 0.1663143. 


X 

Dg ( X,Qo ) 

DZ'(x 

3 Go) 

Df{x 

, Go) 

1 . 0 - 

Iq-02 

1.4844 • 

10+ 00 

1.4650 • 

10 +oo 

6.5605 • 

10 -oi 

5.0- 

Iq-02 

1.1209 ■ 

10 +oo 

1.2757 • 

10 +oo 

6.1705 • 

IO - 01 

1 . 0 - 

IQ-Oi 

8.6823 • 

IQ-Oi 

1.0045 • 

10 +oo 

5.0518- 

10 -oi 

2 . 0 - 

10 -oi 

5.2896 • 

10 -oi 

6.5723 • 

IO - 01 

3.2305 • 

IO - 01 

3.0- 

IO - 01 

2.5532 • 

10 -oi 

4.4560 • 

IO - 01 

1.7382 • 

IO - 01 

4.0- 

IO - 01 

9.4016 - 

Iq-02 

3.0637 • 

IO - 01 

7.7023 • 

Iq-02 

5.0- 

IO - 01 

2.6112 • 

Iq-02 

2.0682 • 

IO - 01 

2.7204 ■ 

Iq-02 

6 . 0 - 

IO ” 01 

5.5517- 

Iq-OS 

1.3142 • 

IO - 01 

7.1061 - 

Iq-OS 

7.0- 

IQ - 01 

1.0569 - 

IQ-OS 

7.4243 • 

l 0 -° 2 

1.1750- 

IQ-OS 

8 . 0 - 

io - 01 

2.1577- 

10 04 

3.3477 • 

l 0 -° 2 

8.5823 • 

10 -os 

9.0- 

IQ-Oi 

2.5554 • 

IQ-OS 

8.6589 • 

IQ-OS 

8.6990 - 

10-07 


Table 6: Same as table 1: FFNS, fi 2 F = 10 2 GeV 2 , jj? R /^ j? f = 2, a s {fi 2 F ) = 0.1525814. 


- 14 - 



































X 

Dg + (x,Q 0 ) 

Df{x 

3 Go) 

Df(x 

,Go) 

1.0- 

10-°2 

2.5036 ■ 

10+oo 

1.8809 ■ 

10+oo 

1.0840 ■ 

10+0° 

5.0- 

Iq-02 

1.6021 ■ 

lO+oo 

1.5851 • 

10+oo 

9.7056 • 

IO- 01 

1.0- 

nr 01 

1.0000 ■ 

lO+oo 

1.1182 • 

10+oo 

6.6513- 

IO- 01 

2.0- 

to- 01 

4.2870 ■ 

10-oi 

6.2699 • 

IO- 01 

3.2715 • 

IO- 01 

3.0- 

to- 01 

1.5797 • 

io- 01 

3.8269 • 

IO- 01 

1.4614 • 

IO- 01 

4.0- 

10-oi 

4.7042 • 

l 0 -° 2 

2.4375 • 

IO- 01 

5.6016- 

Iq-02 

5.0- 

to- 01 

1.1332 • 

l 0 -° 2 

1.5400 ■ 

IO- 01 

1.7502 - 

Iq-02 

6.0- 

10-oi 

2.4107 • 

1Q-0S 

9.1246 • 

l 0 -° 2 

4.0827 - 

l 0 -° 3 

7.0- 

to- 01 

5.4799 • 

10“ 04 

4.7376 • 

l 0 -° 2 

6.0057 - 

IO” 04 

8.0- 

to- 01 

1.2080 ■ 

IO" 04 

1.9061 • 

l 0 -° 2 

3.8052 • 

IQ-OS 

9.0- 

to- 01 

1.2251 ■ 

10-° 5 

4.0856 • 

CO 

o 

1 

O 

2.9598 • 

10-07 


Table 7: Same as table 1: FFNS, fi F = 10 3 GeV 2 , [i 2 R ! [i? F = 1/2, a s ([i 2 F ) = 0.1376668. 


X 

Dg + (x, Qo) 

Df{x 

3 Qo) 

Df{x 

, Go) 

1 . 0 - 

Iq-02 

2.3917- 

10+ 00 

1.9100 • 

lO+oo 

1.1136 • 

lO+oo 

5.0- 

Iq-02 

1.3868 ■ 

lO+oo 

1.4381 • 

10 +oo 

8.1858- 

IO- 01 

1 . 0 - 

IO- 01 

8.9640 ■ 

10 -oi 

1.0374 • 

10 +oo 

5.7848 • 

IO- 01 

2 . 0 - 

IO- 01 

4.1235 • 

10 -oi 

6.1053- 

IO- 01 

3.0621 ■ 

IO- 01 

3.0- 

IO- 01 

1.6110 ■ 

10 -oi 

3.8654 • 

IO- 01 

1.4496 • 

IO- 01 

4.0- 

IQ-Oi 

5.0482 • 

Iq-02 

2.5217- 

10 -oi 

5.8119- 

Iq-02 

5.0- 

IO- 01 

1.2699 • 

Iq-02 

1.6199- 

IO- 01 

1.8821 - 

Iq-02 

6 . 0 - 

IO” 01 

2.7690 ■ 

Iq-OS 

9.7429 • 

l 0 -° 2 

4.5254 • 

Iq-OS 

7.0- 

IQ- 01 

6.2571 - 

10 04 

5.1475 • 

l 0 -° 2 

6.8541 - 

l 0 -° 4 

8 . 0 - 

io- 01 

1.3653 • 

10 04 

2.1204 - 

Iq-02 

4.5049 • 

IQ-OS 

9.0- 

IQ- 01 

1.3932 • 

10 -os 

4.7266 • 

CO 

o 

1 

O 

3.8945 • 

Iq-07 


Table 8: Same as table 1: FFNS, [i 2 F = 10 3 GeV 2 , fj? R /pL 2 F = 1, a s {/j? F ) = 0.1282435. 


X 

Dg ( X,Qo ) 

DZ'(x 

3 Go) 

Df{x 

, Go) 

1 . 0 - 

Iq-02 

2.2362 • 

lO+oo 

1.8597 • 

10 +oo 

1.0617- 

lO+oo 

5.0- 

Iq-02 

1.2724 • 

10 +oo 

1.3614- 

10 +oo 

7.3454 • 

IO - 01 

1 . 0 - 

IQ-Oi 

8.5344 • 

IQ-Oi 

1.0021 • 

10 +oo 

5.3545 • 

10 -oi 

2 . 0 - 

10 -oi 

4.1952 • 

10 -oi 

6.1030- 

IO - 01 

3.0014- 

IO - 01 

3.0- 

IO - 01 

1.7273 • 

10 -oi 

3.9591 - 

IO - 01 

1.4819 - 

IO - 01 

4.0- 

IO - 01 

5.6467 - 

Iq-02 

2.6272 • 

IO - 01 

6.1368- 

Iq-02 

5.0- 

IO - 01 

1.4636 • 

Iq-02 

1.7110- 

IO - 01 

2.0409 ■ 

Iq-02 

6 . 0 - 

IO” 01 

3.2091 ■ 

Iq-OS 

1.0434 • 

IO - 01 

5.0254 • 

Iq-OS 

7.0- 

IQ - 01 

7.1001 - 

10 04 

5.6041 - 

l 0 -° 2 

7.7970 ■ 

l 0 -° 4 

8 . 0 - 

io- 01 

1.5362 • 

10 04 

2.3605 • 

Iq-02 

5.2800 - 

IQ-05 

9.0- 

IQ-Oi 

1.5949 • 

10 -os 

5.4590 - 

Iq-OS 

4.8763 • 

IQ-07 


Table 9: Same as table 1: FFNS, fi 2 F = 10 3 GeV 2 , fi 2 R /fi 2 F = 2, a s (fj, 2 F ) = 0.1200632. 
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X 

Dg + (x,Q 0 ) 

Df{x 

3 Qo ) 

Df(x 

, Qo ) 

1.0- 

10-°2 

3.3009 ■ 

10+oo 

2.3350 • 

10+oo 

1.5488 • 

10+0° 

5.0- 

Iq-02 

1.6386 ■ 

lO+oo 

1.6340 • 

10+oo 

1.0451 ■ 

10+0° 

1.0- 

nr 01 

9.2041 ■ 

10-oi 

1.0909 ■ 

10+oo 

6.6264 • 

IO- 01 

2.0- 

to- 01 

3.3987 • 

10-oi 

5.8002 • 

IO- 01 

2.9899 • 

IO- 01 

3.0- 

to- 01 

1.1232 • 

io- 01 

3.4342 • 

IO- 01 

1.2597- 

IO- 01 

4.0- 

10-oi 

3.1092 • 

l 0 -° 2 

2.1320 ■ 

IO- 01 

4.6024 • 

Iq-02 

5.0- 

to- 01 

7.3648 • 

1Q-OS 

1.3097 • 

IO- 01 

1.3757 - 

Iq-02 

6.0- 

10-oi 

1.6806 • 

1Q-0S 

7.5057 • 

l 0 -° 2 

3.0681 - 

l 0 -° 3 

7.0- 

to- 01 

4.1621 • 

10“ 04 

3.7400 ■ 

IQ- 02 

4.2947 - 

IO” 04 

8.0- 

to- 01 

8.9420 ■ 

10-os 

1.4237 • 

IQ- 02 

2.5637 - 

IQ-OS 

9.0- 

to- 01 

7.9956 • 

l 0 -°6 

2.7874 • 

CO 

o 

1 

O 

1.8560 - 

10-07 


Table 10: Same as table 1: FFNS, fj? F = 10 4 GeV 2 , fj? R /n 2 F = 1/2, a s {^ 2 F ) = 0.1107665. 


X 

Dg + (x,Qo) 

Df{x 

3 Qo) 

Df(x 

,Qo) 

1.0- 

Iq-02 

3.0377 • 

10+ 00 

2.2835 • 

10+oo 

1.4974 • 

lO+oo 

5.0- 

Iq-02 

1.4327 ■ 

lO+oo 

1.4831 • 

10+oo 

8.8904 ■ 

IO- 01 

1.0- 

10- Ql 

8.3511 ■ 

10-oi 

1.0157- 

10+oo 

5.8162 • 

IO- 01 

2.0- 

10- Ql 

3.3048 • 

10-oi 

5.6746 • 

IO- 01 

2.8183- 

IO- 01 

3.0- 

10- Ql 

1.1566- 

10-oi 

3.4830 • 

IO- 01 

1.2567- 

IO- 01 

4.0- 

IQ-Oi 

3.3622 - 

Iq-02 

2.2133 • 

10-oi 

4.7988 • 

Iq-02 

5.0- 

IO- 01 

8.2715 • 

Iq-OS 

1.3823 • 

IO- 01 

1.4861 - 

Iq-02 

6.0- 

IQ-Oi 

1.9166- 

IQ-OS 

8.0443 • 

l 0 -° 2 

3.4158 • 

Iq-OS 

7.0- 

IQ- 01 

4.7207 • 

10 04 

4.0809 - 

l 0 -° 2 

4.9242 • 

l 0 -° 4 

8.0- 

IQ-Oi 

1.0146 • 

10 04 

1.5917- 

l 0 -° 2 

3.0513- 

10-° 5 

9.0- 

IQ- 01 

9.2326 • 

Iq-06 

3.2443 • 

IQ-OS 

2.4600 ■ 

Iq-07 


Table 11: Same as table 1: FFNS, fi 2 F = 10 4 GeV 2 , [i 2 R //j, 2 f = 1, a s (fi 2 F ) = 0.1046608. 


X 

Dg ( X,Qo ) 

DZ'(x 

3 Qo) 

Df{x 

,Qo) 

1.0- 

Iq-02 

2.8153 • 

lO+oo 

2.1928 • 

lO+oo 

1.4047 • 

lO+oo 

5.0- 

Iq-02 

1.3296 • 

10+oo 

1.4068 • 

lO+oo 

8.0505 • 

IO- 01 

1.0- 

IQ-Oi 

8.0483 - 

IQ-Oi 

9.8487 • 

IO- 01 

5.4299 • 

10-oi 

2.0- 

10- Ql 

3.3981 ■ 

10-oi 

5.6968 • 

IO- 01 

2.7806 ■ 

IO- 01 

3.0- 

10- Ql 

1.2511 ■ 

10-oi 

3.5815 • 

IO- 01 

1.2919 - 

IO- 01 

4.0- 

10- Ql 

3.7830 - 

Iq-02 

2.3146 • 

IO- 01 

5.0935 • 

Iq-02 

5.0- 

10- Ql 

9.5276 • 

IQ-OS 

1.4659 • 

IO" 01 

1.6198 • 

Iq-02 

6.0- 

IQ-Oi 

2.2022 • 

IQ-OS 

8.6536 • 

Iq-02 

3.8135 • 

Iq-OS 

7.0- 

IQ- 01 

5.3375 • 

10 04 

4.4661 • 

l 0 -° 2 

5.6342 - 

l 0 -° 4 

8.0- 

IQ-Oi 

1.1502 • 

10 04 

1.7827 • 

Iq-02 

3.6000 ■ 

10-os 

9.0- 

IQ- 01 

1.0725 • 

IQ-OS 

3.7752 • 

Iq-OS 

3.1055 • 

10-07 


Table 12: Same as table 1: FFNS, fi F = IO 4 GeV 2 , ^i 2 R /fj? F = 2, a s ([i 2 F ) = 0.0992090. 
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X 

Dg + (x,Q 0 ) 

Dl + ( x,Q 0 ) 

Ds + ( x,Q 0 ) 

Dc + (x, Qo) 

1.0 • 

10 02 

1.5492 • 

10-cn 

7.6121 • 

10-oi 

-6.0682 

• IO -02 

5.7231 • 

IQ-0! 

5.0 • 

10 02 

9.9725 ■ 

io- 01 

1.2571 - 

10+0° 

5.6325 • 

IO" 01 

4.0424 • 

IO- 01 

1.0 • 

10- 01 

9.4191 ■ 

io- 01 

1.0792 - 

10+0° 

5.4321 • 

IO" 01 

2.2772 ■ 

IO- 01 

2.0 • 

10“ 01 

7.1834- 

IO" 01 

7.4315 • 

io- 01 

3.8246 • 

IO" 01 

7.7980 • 

IO -02 

3.0 • 

10" 01 

4.0469 ■ 

10-oi 

5.1287- 

IO" 01 

2.1550 • 

io- 01 

2.3808 • 

IO -02 

4.0 ■ 

IQ- 01 

1.6712 ■ 

10-oi 

3.5872 • 

io- 01 

9.9051 • 

10-02 

5.7948 ■ 

10-os 

5.0 • 

10~ 01 

4.9921 • 

IO -02 

2.4851 • 

io- 01 

3.6305 • 

Iq-02 

1.0162 • 

10-os 

6.0 • 

10~ 01 

1.0434 • 

l 0 -° 2 

1.6359 • 

io- 01 

9.8915 • 

Iq-OS 

1.0678 • 

IO -04 

7.0 • 

10~ 01 

1.6197- 

10-os 

9.6718 • 

Iq-02 

1.7194 • 

1Q-0S 

3.1017- 

Iq-06 

8.0 • 

10~ 01 

2.8162 • 

IO -04 

4.6376 • 

IQ-02 

1.3318 • 

IO -04 

-5.5107 

• io- 07 

9.0 ■ 

10“ 01 

4.2458 ■ 

10-05 

1.3267 • 

IO -02 

1.3727 • 

10-06 

-5.6280 

• IO -08 


Table 13: Numerical values of gluon, up, strange and charm FFs evolved at NLO in the 
VFNS up to fip = 10 GeV 2 , Mr/m f = 1/2, a s (n F ) = 0.2823897. 


X 

Dg + (x, Qo) 

Dl + ( x,Q 0 ) 

D* + (x,Q 0 ) 

Dc + {x,Qo) 

1.0 • 

l 0 -° 2 

4.9673 ■ 

IO' 01 

9.7661 ■ 

IO" 01 

1.5592 • 

10-oi 

5.4015 ■ 

IQ-O! 

5.0 • 

l 0 -° 2 

8.5963 • 

io- 01 

1.1707- 

10+0° 

4.7431 • 

IO- 01 

3.2641 - 

IO- 01 

1.0 • 

1 0 - 01 

8.2220 • 

1 0 - 01 

1.0056 • 

10+0° 

4.6563 • 

IO- 01 

1.8456 ■ 

IO- 01 

2.0 • 

10~ 01 

6.6765 • 

IO- 01 

7.1528 • 

IO- 01 

3.5174 • 

IO- 01 

6.5025 - 

Iq-02 

3.0 • 

IQ-Oi 

3.9549 • 

IO- 01 

5.0920 • 

IO" 01 

2.0901 • 

IO- 01 

2.0567 - 

Iq-02 

4.0 • 

IQ-Oi 

1.7037 • 

IQ-Oi 

3.6397 ■ 

IO- 01 

9.9961 • 

IO -02 

5.2242 • 

IQ-OS 

5.0 • 

IQ-Oi 

5.2930 • 

l 0 -° 2 

2.5570 ■ 

IO- 01 

3.7780 • 

IO -02 

9.7238 - 

10“ 04 

6.0 • 

IQ-Oi 

1.1495 ■ 

l 0 -° 2 

1.7018 • 

IO" 01 

1.0552 • 

IO -02 

1.1505 ■ 

10“ 04 

7.0 • 

IQ-Oi 

1.8275 ■ 

l 0 -° 3 

1.0181 ■ 

IO- 01 

1.8759 ■ 

10-os 

6.2394 ■ 

10-06 

8.0 • 

IQ-Oi 

3.0556 ■ 

10 04 

4.9605 ■ 

10 o 2 

1.4927 • 

10“ 04 

-9.4928 

■ 10“° 8 

9.0 • 

IQ-Oi 

4.3179 • 

io- 05 

1.4576 • 

10 o 2 

1.6624 • 

Iq-06 

-2.8984 

• IO -08 


Table 14: Same as table 13: VFNS, (j 2 f = 10 GeV 2 , fi 2 R /fi 2 F = 1, a s (^ 2 F ) = 0.2462915. 


X 

Df(x 

,Qo) 

Dl + (x, Qo) 

Ds + ( x,Qo) 

Df {x, Qo) 

1.0 • 

l 0 -° 2 

5.7877 • 

IO- 01 

1.0392 • 

10+0° 

2.1814 • 

io- 01 

4.8959 - 

10 o 1 

5.0 • 

IO" 02 

7.6567 • 

10“ 01 

1.1185 • 

10+0° 

4.1713 ■ 

IO- 01 

2.7920 ■ 

IO" 01 

1.0 • 

IO” 01 

7.5334 • 

1 0 - Ql 

9.6745 • 

IO- 01 

4.2134 • 

IO" 01 

1.5958 - 

IO- 01 

2.0 • 

1 0 - 01 

6.5368 • 

10“ 01 

7.0543 • 

IO- 01 

3.3752 • 

IO- 01 

5.7924 • 

IO -02 

3.0 • 

IO’ 01 

4.0666 ■ 

1 0 - 01 

5.1309 ■ 

IO- 01 

2.0876 ■ 

IO- 01 

1.8923 ■ 

Iq-02 

4.0 • 

IQ-Oi 

1.8218 • 

IQ-Oi 

3.7249 ■ 

IO" 01 

1.0282 • 

IO- 01 

4.9812 - 

10-os 

5.0 • 

1 0 - 01 

5.8534 • 

10 02 

2.6468 • 

IO" 01 

3.9775 • 

IQ-02 

9.6878 - 

IO- 0 4 

6.0 • 

10~ 01 

1.3048 • 

io- 02 

1.7797 • 

IO- 01 

1.1331 • 

Iq-02 

1.2288 • 

IO- 0 4 

7.0 • 

IQ-Oi 

2.0712 • 

CO 

o 

1 

O 

1.0773 • 

IO- 01 

2.0532 • 

IQ-OS 

8.1636 ■ 

10-06 

8.0 • 

IQ-Oi 

3.2831 • 

i 0 -° 4 

5.3326 ■ 

Iq-02 

1.6712 • 

10 - 0 4 

1.6326 ■ 

10-07 

9.0 • 

IQ-Oi 

4.4872 • 

i 0 -° 5 

1.6088 • 

IO -02 

1.9574 • 

Iq-06 

-1.2078 

■ IO -08 


Table 15: Same as table 13: VFNS, /j' 2 f = 10 GeV 2 , [i 2 R /[i 2 F = 2, u s (n 2 F ) = 0.2179093. 


- 17 - 









































X 

Df(x 

,Qo) 

Dl + ( x,Q 0 ) 

Ds + (x, Qo) 

Dl + (x, Qo) 

1.0 • 

l 0 -°2 

1.4392 - 

io+°° 

1.3449 • 

10+°° 

5.3639 ■ 

10-oi 

1.0538 • 

l 0 +oo 

5.0 • 

1q-°2 

1.3984 - 

io+°° 

1.4810 • 

10+°° 

8.3370 • 

10-oi 

6.4282 • 

IO- 01 

1.0 • 

10“ 01 

1.0075 ■ 

io+°° 

1.1281 • 

io+°° 

6.4203 ■ 

IO" 01 

3.3614 • 

IO" 01 

2.0 • 

10~ 01 

5.3182 • 

10~ 01 

6.8312 ■ 

IO" 01 

3.5870 • 

10-oi 

1.0195 ■ 

IO- 01 

3.0 • 

10~ 01 

2.2943 • 

io- 01 

4.3596 • 

IO" 01 

1.7444 • 

IO" 01 

2.7940 • 

Iq-02 

4.0 • 

10“ 01 

7.7030 • 

Iq-02 

2.8767 • 

10-oi 

7.1576 ■ 

l 0 -°2 

6.1888 • 

IQ-OS 

5.0 • 

10~ 01 

1.9870 • 

Iq-02 

1.8857 • 

IO" 01 

2.3781 • 

l 0 -°2 

1.0124 • 

IQ-OS 

6.0 • 

10~ 01 

4.0894 • 

CO 

o 

1 

O 

1.1666 • 

10-oi 

5.8981 • 

IQ-OS 

1.0872 • 

IO -04 

7.0 • 

10” 01 

8.2563 • 

10 04 

6.3903 • 

l 0 -°2 

9.2798 • 

10“ 04 

6.8820 • 

10-06 

8.0 • 

10” 01 

1.8409 ■ 

IQ- 04 

2.7637 ■ 

l 0 -°2 

6.3832 • 

10 -os 

3.3464 • 

10-07 

9.0 • 

10” 01 

2.2470 • 

10~ 05 

6.6713 ■ 

1Q-0S 

5.5920 • 

Iq-07 

-3.5677 

. 10-10 


Table 16: Same as table 13: VFNS, fi 2 F = 10 2 GeV 2 , jj? R /fx 2 F = 1/2, a s (y? F ) = 0.1938346. 


X 

Df(x 

,Qo) 

Dl + (x,Q 0 ) 

Ds + (x, Qo) 

Dc + ( x,Qo ) 

1.0 • 

Iq-02 

1.5355 • 

10 +0° 

1.4712 • 

io+°° 

6.6365 • 

IO" 01 

9.7376 ■ 

IQ-O! 

5.0 • 

Iq-02 

1.2070 • 

10 +0° 

1.3548 • 

io+°° 

7.0408 • 

IO - 01 

5.3020 ■ 

IO" 01 

1.0 • 

IQ-Oi 

8.9152 • 

IQ-Oi 

1.0465 ■ 

io+°° 

5.5597 ■ 

IO - 01 

2.8080 ■ 

IO" 01 

2.0 • 

IQ-Oi 

5.0053 • 

IQ-Oi 

6.6057 • 

IO - 01 

3.3280 ■ 

IO - 01 

8.8346 ■ 

Iq-02 

3.0 • 

IQ-Oi 

2.2730 ■ 

IQ-Oi 

4.3605 ■ 

IO" 01 

1.7082 ■ 

IO - 01 

2.5184 • 

Iq-02 

4.0 • 

IQ-Oi 

7.9809 ■ 

Iq-02 

2.9423 ■ 

IO - 01 

7.3049 ■ 

Iq-02 

5.8344 • 

IQ-OS 

5.0 • 

IQ-Oi 

2.1440 • 

Iq-02 

1.9573 • 

IO - 01 

2.5068 • 

Iq-02 

1.0131 • 

IQ-OS 

6.0 • 

IQ-Oi 

4.5498 • 

Iq-OS 

1.2260 • 

IO - 01 

6.3849 • 

Iq-OS 

1.2059 • 

IO -04 

7.0 • 

IQ-Oi 

9.1769 • 

1Q 04 

6.8098 • 

IQ-02 

1.0297 ■ 

IQ-OS 

9.6284 • 

10-06 

8.0 • 

IQ-Oi 

1.9945 • 

1 0 04 

3.0013 • 

Iq-02 

7.3001 • 

10-os 

6.6294 • 

10-07 

9.0 ■ 

IQ-Oi 

2.3977 • 

IQ-OS 

7.4764 • 

Iq-OS 

6.9708 • 

10-07 

1.5982 • 

Iq-08 


Table 17: Same as table 13: VFNS, n 2 F = 10 2 GeV 2 , n 2 R /n 2 F = 1, a s (/4) = 0.1777432. 


X 

Df{x 

,Qo) 

DZ + (x,Q 0 ) 

Ds + (x,Qo) 

Dc + {x,Qo) 

1.0 • 

Iq-02 

1.4786 • 

10+0° 

1.4728 • 

io+°° 

6.6434 • 

IO" 01 

8.8195 - 

IQ-O! 

5.0 • 

Iq-02 

1.0919 • 

10+0° 

1.2820 • 

10+°° 

6.2492 • 

IO- 01 

4.6372 • 

IO- 01 

1.0 • 

IQ-Oi 

8.3428 • 

IQ-Oi 

1.0066 ■ 

io+°° 

5.0888 • 

IO- 01 

2.5032 • 

10-oi 

2.0 • 

IQ-Oi 

5.0011 ■ 

IQ-Oi 

6.5597 • 

IO- 01 

3.2292 ■ 

IO- 01 

8.1847- 

Iq-02 

3.0 • 

IQ-Oi 

2.3909 ■ 

IQ-Oi 

4.4351 • 

IO- 01 

1.7285 ■ 

IO- 01 

2.4214 • 

iq-02 

4.0 • 

IQ-Oi 

8.7547 • 

Iq-02 

3.0431 • 

IO- 01 

7.6299 ■ 

IQ-02 

5.8296 ■ 

IQ-OS 

5.0 • 

IQ-Oi 

2.4307 ■ 

Iq-02 

2.0505 ■ 

IO- 01 

2.6862 • 

Iq-02 

1.0577 - 

IQ-OS 

6.0 • 

IQ-Oi 

5.2310 • 

Iq-OS 

1.3004 • 

IO- 01 

6.9967 ■ 

Iq-OS 

1.3345 - 

10“ 04 

7.0 • 

IQ-Oi 

1.0313 • 

IQ-OS 

7.3297 • 

IQ-02 

1.1537- 

IQ-OS 

1.1638 • 

10-os 

8.0 • 

IQ-Oi 

2.1805 • 

10 04 

3.2950 • 

Iq-02 

8.4013 • 

10-os 

8.8603 - 

10-07 

9.0 • 

IQ-Oi 

2.6290 • 

IQ-OS 

8.4809 • 

Iq-OS 

8.5046 • 

10-07 

2.8561 - 

IQ-OS 


Table 18: Same as table 13: VFNS, ji F = 10 2 GeV 2 , [i 2 R /n 2 F = 2, ot s (y? F ) = 0.1637695. 
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X 

Df(x 

,Qo) 

DZ + ( x,Q 0 ) 

DZ + ( x,Q 0 ) 

DZ + {x, Qo) 

1.0 • 

l 0 -° 2 

2.5055 • 

10 +°° 

1.9150 ■ 

io+°° 

1.1192 • 

10 +°° 

1.5585 • 

lO+oo 

5.0 • 

10” 02 

1.5215 • 

10 +°° 

1.5863 ■ 

io+°° 

9.7389 ■ 

IO" 01 

7.6647 • 

IO" 01 

1.0 • 

10” 01 

9.3252 • 

IO" 01 

1.1145 • 

io+°° 

6.6353 • 

IO" 01 

3.6979 ■ 

io- 01 

2.0 • 

10” 01 

3.9236 ■ 

10” 01 

6.2284 • 

io- 01 

3.2456 ■ 

IO" 01 

1.0074 • 

io- 01 

3.0 • 

10~ 01 

1.4328 • 

10~ 01 

3.7950 • 

io- 01 

1.4449 ■ 

IO" 01 

2.5391 • 

Iq-02 

4.0 • 

10~ 01 

4.2628 • 

IQ-02 

2.4134 • 

IO" 01 

5.5238 • 

10-02 

5.2514 • 

10 -os 

5.0 • 

10~ 01 

1.0403 • 

1 Q-Q2 

1.5220 • 

IO" 01 

1.7218 • 

Iq-02 

8.1976 ■ 

10-04 

6.0 • 

10~ 01 

2.2992 • 

l 0 -° 3 

8.9983 • 

l 0 -°2 

4.0074 • 

IO -03 

8.9814 • 

10-05 

7.0 • 

10" 01 

5.4993 • 

10” 04 

4.6600 • 

IQ-02 

5.8831 • 

10-04 

7.3959 • 

10-06 

8.0 • 

10- 01 

1.2453 ■ 

10 04 

1.8685 • 

Iq-02 

3.7254 ■ 

10-05 

6.0313 ■ 

10-07 

9.0 ■ 

10” 01 

1.2788 • 

IO” 05 

3.9843 • 

IO -03 

2.9477 • 

IQ-07 

1.4005 ■ 

10-08 


Table 19: Same as table 13: VFNS, fi 2 F = 10 3 GeV 2 , jj? R /fx 2 F = 1/2, a s (y? F ) = 0.1501564. 


X 

Df{x 

,Qo) 

DZ + {x,Q 0 ) 

DZ + (z,Qo) 

DZ + {x, Qo) 

1.0 • 

Iq-02 

2.3873 • 

10 +0° 

1.9405 • 

io+°° 

1.1453- 

io+°° 

1.4044 - 

lO+oo 

5.0 • 

Iq-02 

1.3242 • 

10 +0° 

1.4453 ■ 

io+°° 

8.2905 ■ 

IO- 01 

6.4117- 

10 -oi 

1.0 • 

IQ-Oi 

8.3590 ■ 

IQ-Oi 

1.0361 • 

io+°° 

5.8054 • 

IO- 01 

3.1442 • 

10-04 

2.0 • 

IQ-Oi 

3.7371 ■ 

IQ-Oi 

6.0550 • 

IO- 01 

3.0363 ■ 

IO- 01 

8.9037 • 

Iq-02 

3.0 • 

IQ-Oi 

1.4356 ■ 

IQ-Oi 

3.8162 • 

IO- 01 

1.4253 • 

IO- 01 

2.3359 ■ 

Iq-02 

4.0 • 

IQ-Oi 

4.4636 ■ 

Iq-02 

2.4804 • 

10 -oi 

5.6760 ■ 

Iq-02 

5.0552 - 

IQ-OS 

5.0 • 

IQ-Oi 

1.1303 • 

Iq-02 

1.5874 • 

IO- 01 

1.8272 • 

Iq-02 

8.3751 - 

Iq-04 

6.0 • 

IQ-Oi 

2.5472 • 

1Q 03 

9.5071 • 

Iq-02 

4.3682 • 

Iq-OS 

1.0082 • 

Iq-04 

7.0 • 

IQ-Oi 

6.0528 • 

1 Q 04 

4.9972 • 

Iq-02 

6.5773 • 

Iq-04 

9.5737 - 

Iq-06 

8.0 • 

IQ-Oi 

1.3583 • 

10 04 

2.0445 • 

Iq-02 

4.2971 • 

IQ-05 

8.4913 - 

Iq-07 

9.0 ■ 

IQ-Oi 

1.4014 ■ 

Iq-05 

4.5074 • 

Iq-OS 

3.7152 • 

Iq-07 

2.5215 ■ 

Iq-08 


Table 20: Same as table 13: VFNS, /jl 2 f = 10 3 GeV 2 , fx R / = 1, a s {^i 2 F ) = 0.1404546. 


X 

Df{x 

,Qo) 

DZ + {x,Q 0 ) 

DZ + (x,Qo) 

DZ + {x,Qo) 

1.0 • 

Iq-02 

2.2233 - 

io+°° 

1.8842 • 

io+°° 

1.0876 • 

10 +0° 

1.2672 • 

lO+oo 

5.0 • 

Iq-02 

1.2150 • 

10 +°° 

1.3672 • 

10 +°° 

7.4390 • 

10-°4 

5.6957 • 

10-04 

1.0 • 

Iq-04 

7.9516 ■ 

Iq-04 

1.0004 • 

io+°° 

5.3729 • 

10-04 

2.8555 • 

10-04 

2.0 • 

Iq-04 

3.7904 ■ 

iq-04 

6.0485 ■ 

10-°4 

2.9735 ■ 

10-04 

8.4231 • 

IQ-02 

3.0 • 

Iq-04 

1.5315 ■ 

iq-04 

3.9051 • 

10-°4 

1.4548 ■ 

10-04 

2.2956 ■ 

Iq-02 

4.0 • 

iq-04 

4.9562 • 

iq-02 

2.5810 • 

10-04 

5.9804 • 

IQ-02 

5.1639 ■ 

IQ-OS 

5.0 • 

iq-04 

1.2888 • 

Iq-02 

1.6741 • 

10-04 

1.9758 • 

iq-02 

8.9297 • 

Iq-04 

6.0 • 

iq-04 

2.9074 • 

1Q 03 

1.0162 • 

10-04 

4.8342 • 

IQ-OS 

1.1312 • 

Iq-04 

7.0 • 

iq-04 

6.7562 • 

1Q 04 

5.4272 • 

IQ-02 

7.4501 • 

Iq-04 

1.1337- 

10 -os 

8.0 • 

iq-04 

1.5053 • 

1Q 04 

2.2688 • 

iq-02 

5.0079 • 

10 -os 

1.0444 • 

10-06 

9.0 • 

iq-04 

1.5798 • 

Iq-05 

5.1832 • 

IQ-OS 

4.6023 • 

10-07 

3.5222 • 

Iq-08 


Table 21: Same as table 13: VFNS, /i F = 10 3 GeV 2 , [i 2 R /n 2 F = 2, a s (y? F ) = 0.1316906. 
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X 

Df(x 

,Qo) 

DZ + {x,Q 0 ) 

DZ + ( x,Q 0 ) 

DZ + {x, Qo) 

1.0 • 

l 0 -° 2 

3.2870 • 

10 +°° 

2.3905 ■ 

io+°° 

1.6064 • 

10 +0° 

1.9877 • 

lO+oo 

5.0 • 

10~ 02 

1.5277 • 

10 +°° 

1.6309 ■ 

io+°° 

1.0465 • 

10 +0° 

8.2966 ■ 

IO" 01 

1.0 • 

10- 01 

8.3390 • 

10“ 01 

1.0812 • 

io+°° 

6.5709 ■ 

IO- 01 

3.7534 • 

IO" 01 

2.0 • 

10” 01 

2.9811 ■ 

10” 01 

5.7123 ■ 

io- 01 

2.9337 ■ 

IO- 01 

9.4061 • 

Iq-02 

3.0 • 

10~ 01 

9.6880 • 

1 Q-Q2 

3.3697 • 

io- 01 

1.2270 • 

IO- 01 

2.2257 • 

Iq-02 

4.0 • 

10~ 01 

2.6760 • 

l 0 -°2 

2.0844 • 

IO" 01 

4.4550 • 

Iq-02 

4.3784 • 

Iq-OS 

5.0 • 

10~ 01 

6.4771 • 

IQ-OS 

1.2750 • 

IO" 01 

1.3238 • 

Iq-02 

6.6282 • 

IO -04 

6.0 • 

10~ 01 

1.5557 ■ 

IQ-OS 

7.2690 • 

l 0 -°2 

2.9347 ■ 

Iq-OS 

7.4275 • 

10-05 

7.0 • 

10" 01 

4.0416 ■ 

10~ 04 

3.5993 • 

IQ -02 

4.0828 • 

l 0 -° 4 

7.0423 • 

IQ -06 

8.0 • 

10- 01 

8.8370 ■ 

IQ-OS 

1.3588 • 

Iq-02 

2.4246 • 

10-05 

6.4291 - 

Iq-07 

9.0 ■ 

10” 01 

7.9382 ■ 

Iq-06 

2.6253 ■ 

Iq-OS 

1.7799 ■ 

IQ -07 

1.6186 ■ 

Iq-08 


Table 22: Same as table 13: VFNS, [i 2 F = 10 4 GeV 2 , fi 2 R /n 2 F = 1/2, a s (fx F ) = 0.1228358. 


X 

Df(x 

,Qo) 

DZ + {x,Qo) 

DZ + (z,Qo) 

DZ + {x, Qo) 

1.0 • 

Iq-02 

3.0177- 

10 +0° 

2.3322 • 

io+°° 

1.5485 • 

io+°° 

1.7693 • 

lO+oo 

5.0 • 

Iq-02 

1.3404 • 

10 +0° 

1.4855 ■ 

io+°° 

8.9706 ■ 

IO- 01 

7.0101 • 

IO" 01 

1.0 • 

IQ-Oi 

7.5535 ■ 

IQ-Oi 

1.0082 - 

io+°° 

5.7953 ■ 

IO- 01 

3.2288 - 

IO" 01 

2.0 • 

IQ-Oi 

2.8672 ■ 

IQ-Oi 

5.5778 • 

IO- 01 

2.7622 ■ 

IO- 01 

8.4158 - 

Iq-02 

3.0 • 

IQ-Oi 

9.7934 • 

Iq-02 

3.4018 • 

IO- 01 

1.2168 • 

IO- 01 

2.0730 ■ 

Iq-02 

4.0 • 

IQ-Oi 

2.8222 • 

Iq-02 

2.1494 ■ 

10 -oi 

4.5993 ■ 

Iq-02 

4.2677 - 

IQ-OS 

5.0 • 

IQ-Oi 

7.0525 ■ 

Iq-OS 

1.3341 • 

IO- 01 

1.4109 • 

Iq-02 

6.8549 - 

Iq-04 

6.0 • 

IQ-Oi 

1.7140 • 

IQ-OS 

7.7074 • 

Iq-02 

3.2124 • 

Iq-OS 

8.3824 • 

IQ-OS 

7.0 • 

IQ-Oi 

4.4339 • 

1 Q 04 

3.8755 • 

Iq-02 

4.5846 • 

Iq-04 

8.8037 - 

Iq-06 

8.0 • 

IQ-Oi 

9.7044 • 

Iq-05 

1.4938 • 

Iq-02 

2.8103 • 

IQ-05 

8.3715 - 

Iq-07 

9.0 ■ 

IQ-Oi 

8.8491 ■ 

Iq-06 

2.9872 • 

Iq-OS 

2.2567 ■ 

Iq-07 

2.4436 ■ 

Iq-08 


Table 23: Same as table 13: VFNS, /jl 2 f = 10 4 GeV 2 , fx R / = 1, a s {^i 2 F ) = 0.1163266. 


X 

Df{x 

,Qo) 

DZ + {x, Qo) 

DZ + (x,Qo) 

DZ + {x,Q 0 ) 

1.0 • 

Iq-02 

2.7851 • 

10 +0° 

2.2317 - 

10 +0° 

1.4460 - 

10 +0° 

1.5968 • 

lO+oo 

5.0 • 

Iq-02 

1.2438 • 

10 +0° 

1.4079 - 

10 +0° 

8.1203 ■ 

IO- 01 

6.2973 - 

IO- 01 

1.0 • 

IQ-Oi 

7.2742 • 

IQ-Oi 

9.7721 • 

IO- 01 

5.4096 ■ 

IO- 01 

2.9692 • 

10 -oi 

2.0 • 

IQ-Oi 

2.9399 ■ 

IQ-Oi 

5.5973 ■ 

IO- 01 

2.7236 ■ 

IO- 01 

8.0665 ■ 

IQ-02 

3.0 • 

IQ-Oi 

1.0543 • 

IQ-Oi 

3.4959 ■ 

IO- 01 

1.2494 • 

IO- 01 

2.0642 - 

iq-02 

4.0 • 

IQ-Oi 

3.1514- 

Iq-02 

2.2459 ■ 

IO- 01 

4.8735 ■ 

Iq-02 

4.4159 ■ 

IQ-OS 

5.0 • 

IQ-Oi 

8.0281 ■ 

IQ-OS 

1.4131 • 

IO- 01 

1.5344 • 

Iq-02 

7.3941 • 

IQ-04 

6.0 • 

IQ-Oi 

1.9391 • 

IQ-OS 

8.2795 • 

Iq-02 

3.5763 ■ 

Iq-OS 

9.4646 - 

10-05 

7.0 • 

IQ-Oi 

4.9426 ■ 

1 Q 04 

4.2336 • 

Iq-02 

5.2269 ■ 

Iq-04 

1.0331 - 

IQ-05 

8.0 • 

IQ-Oi 

1.0864 • 

1 Q 04 

1.6691 • 

Iq-02 

3.2993 • 

IQ-05 

1.0083 - 

Iq-06 

9.0 • 

IQ-Oi 

1.0149 • 

Iq-05 

3.4642 • 

Iq-OS 

2.8196 ■ 

Iq-07 

3.2481 - 

Iq-08 


Table 24: Same as table 13: VFNS, ji F = 10 4 GeV 2 , [i 2 R /n 2 F = 2, a s (y? F ) = 0.1102998. 
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X 

Df(x 

,Qo) 

Df{x 

, Qo) 

D* + (x, Qo) 

1.0- 

to -02 

-1.8792 

- 10 +o ° 

-4.7902 

. io-o 1 

-1.3037 

• 10+°° 

5.0- 

10-02 

6.9280 • 

10-oi 

1.0172 • 

lO+oo 

3.1863 • 

10- 01 

1.0- 

10- 01 

9.1279 ■ 

10-oi 

1.0561 • 

l 0 +oo 

5.1619 ■ 

10- 01 

2.0- 

10-oi 

7.5327 • 

10-oi 

7.7432 • 

10-oi 

4.1056 ■ 

10- 01 

3.0- 

10-oi 

4.2646 ■ 

10-oi 

5.3095 • 

10- Ql 

2.3102 • 

10-oi 

4.0- 

10-oi 

1.7559 ■ 

10-oi 

3.6539 • 

10-oi 

1.0401 • 

10-oi 

5.0- 

10-oi 

5.2453 • 

IQ-02 

2.5022 • 

10-oi 

3.7271 • 

10-02 

6.0- 

10-oi 

1.1152 • 

IQ-02 

1.6363 • 

10-oi 

9.9634 • 

1Q-0S 

7.0- 

10-oi 

1.8531 • 

IQ-OS 

9.6265 • 

Iq-02 

1.7089 • 

1Q-0S 

8.0- 

10-oi 

3.5067 ■ 

1Q- 04 

4.5884 • 

l 0 -°2 

1.3207 • 

10“ 04 

9.0- 

10-oi 

5.2702 • 

IQ-OS 

1.2993 • 

Iq-02 

1.5603 ■ 

10-06 


Table 25: Numerical values of gluon, up and strange FFs evolved at NNLO in the FFNS 
up to ii 2 f = 10 GeV 2 , fi 2 R /n 2 F = 1/2, a s (y? F ) = 0.2749250. 


X 

Df(x 

> Qo) 

Df{x 

,Qo) 

Df(x 

> Qo) 

1.0- 

Iq-02 

-5.2241 

- io- 01 

3.5771 - 

IQ-Oi 

-4.6473 

■ io- 01 

5.0- 

Iq-02 

8.0768 • 

IO- 01 

1.1145 • 

lO+oo 

4.1771 - 

IQ-Oi 

1.0- 

10-oi 

8.7624 • 

IO- 01 

1.0316 ■ 

lO+oo 

4.9234 • 

IQ-Oi 

2.0- 

10- 01 

7.1263 • 

IO- 01 

7.4120 ■ 

IQ-Oi 

3.7806 • 

IQ-Oi 

3.0- 

10-oi 

4.1162 • 

IO- 01 

5.1782 • 

IQ-Oi 

2.1820 ■ 

IQ-Oi 

4.0- 

io- 01 

1.7300 • 

10-oi 

3.6346 • 

IQ-Oi 

1.0125 • 

IQ-Oi 

5.0- 

10- 01 

5.2577 • 

IQ-02 

2.5221 - 

IQ-Oi 

3.7282 • 

Iq-02 

6.0- 

10- 01 

1.1250 • 

iq-02 

1.6637 • 

IQ-Oi 

1.0191 - 

Iq-02 

7.0- 

10- 01 

1.8117- 

IQ-OS 

9.8662 • 

Iq-02 

1.7793 • 

Iq-OS 

8.0- 

10-os 

3.2016- 

o 

1 

o 

4.7515 • 

Iq-02 

1.3923 • 

o 

1 

o 

9.0- 

10- 01 

4.6958 • 

IQ-OS 

1.3691 - 

Iq-02 

1.5652 • 

IQ-06 


Table 26: Same as table 25: FFNS, fj? F = 10 GeV 2 , n 2 R /fi 2 F = 1, a s (^ 2 F ) = 0.2369645. 


X 

Df(x 

, Qo) 

Df{x 

! Qo) 

Df(x 

, Qo) 

1.0- 

Iq-02 

-2.1326 

• IO -02 

6.7241 - 

IQ-Oi 

-1.4945 

■ io- 01 

5.0- 

l 0 -° 2 

7.9209 - 

IQ-Oi 

1.1159 ■ 

1Q+00 

4.1707 • 

IQ-Oi 

1.0- 

10- 01 

8.2704 - 

IQ-Oi 

1.0026 ■ 

lO+oo 

4.6041 - 

IQ-Oi 

2.0- 

10- 01 

6.9265 • 

IQ-Oi 

7.2545 • 

IQ-Oi 

3.6024 • 

IQ-Oi 

3.0- 

10- 01 

4.1308 • 

IQ-Oi 

5.1581 - 

IQ-Oi 

2.1411 - 

IQ-Oi 

4.0- 

10-oi 

1.7841 • 

IQ-Oi 

3.6738 • 

IQ-Oi 

1.0192 • 

io- 01 

5.0- 

10- 01 

5.5504 • 

1Q 02 

2.5756 • 

IQ-Oi 

3.8335 • 

1Q 02 

6.0- 

10-oi 

1.2063 • 

1Q 02 

1.7131 - 

IQ-Oi 

1.0668 • 

1Q 02 

7.0- 

10- 01 

1.9185 • 

CO 

o 

1 

o 

1.0250 • 

IQ-Oi 

1.8924 • 

IQ-OS 

8.0- 

io- 01 

3.2144 • 

o 

1 

O 

4.9950 • 

10 02 

1.5055 • 

O 

1 

O 

9.0- 

10- 01 

4.5968 • 

IQ-OS 

1.4680 • 

I— 1 

o 

1 

o 

to 

1.7142 • 

IQ-oe 


Table 27: Same as table 25: FFNS, /j 2 f = 10 GeV 2 , n 2 R /fi 2 F = 2, a s ([i 2 F ) = 0.2087674. 
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X 

Df(x 

3 Go) 

DZ + {x,Q 0 ) 

Df{x 

> Go) 

1.0 • 

10~ 02 

-9.2616 

■ io-oi 

-2.4841 

. lo-oi 

-1.0610 

. 10+°° 

5.0 • 

Iq-02 

1.2010 • 

io+°° 

1.2584 • 

io+°° 

6.0423 • 

IQ-Oi 

1.0 • 

10“ 01 

1.0385 • 

io+°° 

1.1216 ■ 

io+°° 

6.2957 • 

IQ-Oi 

2.0 • 

10“ 01 

5.8416 ■ 

10~ 01 

7.1622 • 

IQ-Oi 

3.8735 • 

IQ-Oi 

3.0 • 

10~ 01 

2.5511 • 

10~ 01 

4.5378 • 

IQ-Oi 

1.8854 • 

IQ-Oi 

4.0 • 

10~ 01 

8.6015 • 

IQ-02 

2.9513 • 

IQ-Oi 

7.6046 • 

Iq-02 

5.0 • 

10~ 01 

2.2294 • 

iq-02 

1.9170 • 

IQ-Oi 

2.4786 • 

Iq-02 

6.0 • 

10~ 01 

4.6376 ■ 

1Q-0S 

1.1809 • 

IQ-Oi 

6.0531 • 

Iq-OS 

7.0 • 

IQ- 01 

9.4689 • 

1 Q 04 

6.4518 • 

IQ-02 

9.4394 • 

1 Q 04 

8.0 • 

10“ 01 

2.0893 ■ 

1 Q 04 

2.7820 • 

iq-02 

6.5371 ■ 

IQ-OS 

9.0 • 

10” 01 

2.4937 ■ 

1Q-0S 

6.6798 • 

1Q-OS 

6.9517- 

IQ-07 


Table 28: Same as table 25: FFNS, / l 2 f = 10 2 GeV 2 , n 2 R //j' 2 f = 1/2, a s {^? F ) = 0.1808461. 


X 

Df(x 

j Go) 

Df{x 

,Go) 

Df{x 

,Go) 

1.0- 

Iq-02 

4.3018 • 

IQ-Oi 

7.0077 • 

10-oi 

-1.0938 

- io- 01 

5.0- 

Iq-02 

1.2356 • 

lO+oo 

1.3198 • 

1Q+00 

6.6718 • 

IQ-Oi 

1.0- 

10-oi 

9.8493 • 

IQ-Oi 

1.0845 • 

1Q+00 

5.9300 ■ 

IQ-Oi 

2.0- 

10- 01 

5.5526 • 

IQ-Oi 

6.8698 • 

IQ-Oi 

3.5853 • 

IQ-Oi 

3.0- 

10-oi 

2.4764 • 

IQ-Oi 

4.4447 • 

IQ-Oi 

1.7907 • 

IQ-Oi 

4.0- 

10-°! 

8.5205 • 

IQ-02 

2.9474 • 

IQ-Oi 

7.4407 • 

Iq-02 

5.0- 

10- 01 

2.2425 • 

iq-02 

1.9390 ■ 

IQ-Oi 

2.4911 - 

Iq-02 

6.0- 

10- 01 

4.6706 - 

CO 

o 

1 

o 

1.2050 - 

IQ-Oi 

6.2193- 

Iq-OS 

7.0- 

10- 01 

9.3228 • 

1 Q 04 

6.6410 ■ 

Iq-02 

9.8675 • 

1 Q 04 

8.0- 

10- 01 

2.0213 • 

1 Q 04 

2.8963 • 

Iq-02 

6.9032 - 

IQ-OS 

9.0- 

10- 01 

2.4256 • 

IQ-OS 

7.0874 • 

Iq-OS 

6.8263 • 

IQ-07 


Table 29: Same as table 25: FFNS, /j 2 f = 10 2 GeV 2 , Mr/Mf = 1, a s {^ 2 F ) = 0.1644438. 


X 

Dg ( X,Qo ) 

DZ"(x 

3 Go) 

Df{x 

, Go) 

1.0- 

Iq-02 

8.7835 • 

IO- 01 

1.0344 • 

lO+oo 

2.2467 - 

10-oi 

5.0- 

Iq-02 

1.1872 • 

lO+oo 

1.3021 • 

10+0° 

6.4691 - 

IQ-Oi 

1.0- 

10-°! 

9.3764 • 

IO- 01 

1.0515 • 

lQ + °° 

5.5676 • 

IQ-Oi 

2.0- 

IO- 01 

5.4582 • 

io- 01 

6.7498 • 

IQ-Oi 

3.4415 • 

IQ-Oi 

3.0- 

IO- 01 

2.5133- 

10-oi 

4.4508 • 

IQ-Oi 

1.7702 • 

IQ-Oi 

4.0- 

IO- 01 

8.8871 ■ 

Iq-02 

2.9953 • 

IQ-Oi 

7.5476 • 

Iq-02 

5.0- 

IO- 01 

2.3895 • 

Iq-02 

1.9914 • 

IQ-Oi 

2.5824 • 

Iq-02 

6.0- 

IO" 01 

5.0127- 

Iq-OS 

1.2489 • 

IQ-Oi 

6.5675 • 

Iq-OS 

7.0- 

IO- 01 

9.8118- 

IO -04 

6.9528 • 

Iq-02 

1.0596 • 

IQ-OS 

8.0- 

IO" 01 

2.0853 • 

IO -04 

3.0739 • 

Iq-02 

7.5407 - 

IQ-OS 

9.0- 

IO- 01 

2.5117- 

10-° 5 

7.6942 • 

Iq-OS 

7.4918 • 

Iq-07 


Table 30: Same as table 25: FFNS, ^l 2 f = 10 2 GeV 2 , [i 2 R /fi 2 F = 2, a s (fi 2 F ) = 0.1508899. 
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X 

Df(x 

, Go) 

Df(x 

,Go) 

Df{x 

,Go) 

1 .0- 

l 0 -°2 

1.6184 • 

io- 01 

1.8297 • 

IO- 01 

-6.1839 

• IO" 01 

5.0- 

l 0 -°2 

1.4257 • 

lO+oo 

1.3940 ■ 

10 +oo 

7.7217- 

IO- 01 

1 .0- 

to- 01 

1.0082 • 

lO+oo 

1.1231 ■ 

lO+oo 

6.6379 - 

IO- 01 

2 .0- 

to- 01 

4.5375 • 

IO- 01 

6.5868 • 

IO- 01 

3.5427 • 

IO- 01 

3.0- 

to- 01 

1.6934 • 

IO” 01 

3.9875 • 

IO- 01 

1.5840 ■ 

IO- 01 

4.0- 

to- 01 

5.0877 • 

l 0 -° 2 

2.5062 • 

IO- 01 

5.9787 - 

Iq-02 

5.0- 

to- 01 

1.2428 • 

l 0 -° 2 

1.5713- 

IO- 01 

1.8366 • 

Iq-02 

6 .0- 

to- 01 

2.7101 • 

CO 

o 

1 

O 

9.2822 • 

Iq-02 

4.2298 • 

l 0 -° 3 

7.0- 

10 -oi 

6.3006 • 

i 0 -° 4 

4.8139 - 

Iq-02 

6.1899- 

IO” 04 

8 .0- 

10 -oi 

1.3909 • 

i 0 -° 4 

1.9346 • 

Iq-02 

3.9802 - 

IQ-OS 

9.0- 

IQ" 01 

1.4025 • 

10-° 5 

4.1386- 

IO” 03 

3.9353 • 

10-07 


Table 31: Same as table 25: FFNS, fi 2 F = 10 3 GeV 2 , ji 2 R j/ j' f = 1/2, a s {^? F ) = 0.1361858. 


X 

Dg + (x,Qo) 

Df{x 

3 Go) 

Df{x 

, Go) 

1.0- 

Iq-02 

1.3748 • 

lO+oo 

1.1309 • 

10+oo 

3.3200 ■ 

10-oi 

5.0- 

Iq-02 

1.4172 ■ 

lO+oo 

1.4306 ■ 

10+oo 

8.1015 • 

IO - 01 

1.0- 

IO - 01 

9.5459 ■ 

10-oi 

1.0821 • 

10+oo 

6.2332 • 

IO - 01 

2.0- 

10-oi 

4.3327 • 

10-oi 

6.3347 • 

IO - 01 

3.2929 - 

IO - 01 

3.0- 

IO - 01 

1.6513- 

10-oi 

3.9183- 

IO - 01 

1.5105 • 

IO - 01 

4.0- 

IQ-Oi 

5.0591 ■ 

Iq-02 

2.5091 - 

10-oi 

5.8690 ■ 

Iq-02 

5.0- 

IO - 01 

1.2527- 

10 02 

1.5924 • 

IO - 01 

1.8506 - 

Iq-02 

6.0- 

IO” 01 

2.7322 • 

IQ-OS 

9.4895 • 

l 0 -° 2 

4.3543 • 

Iq-OS 

7.0- 

IQ - 01 

6.2833 • 

10 04 

4.9650 • 

l 0 -° 2 

6.4788 • 

l 0 -° 4 

8.0- 

io- 01 

1.3864 • 

10 04 

2.0188- 

l 0 -° 2 

4.2019 • 

10-os 

9.0- 

IQ - 01 

1.4184 • 

10-os 

4.4034 • 

CO 

o 

1 

O 

3.8275 • 

Iq-07 


Table 32: Same as table 25: FFNS, /j 2 f = 10 3 GeV 2 , ji 2 R //j? f = 1, a s {^ 2 F ) = 0.1269012. 


X 

Dg ( X,Qo ) 

DZ"(x 

3 Go) 

Df{x 

, Go) 

1.0- 

Iq-02 

1.7277 • 

10+ot) 

1.4412 • 

lO+oo 

6.4263 • 

10-oi 

5.0- 

Iq-02 

1.3585 • 

10+oo 

1.4033 • 

lO+oo 

7.8017- 

IO- 01 

1.0- 

IQ-Oi 

9.1554- 

IQ-Oi 

1.0502 • 

10+oo 

5.8810 ■ 

10-oi 

2.0- 

10-oi 

4.2926 • 

10-oi 

6.2449 • 

IO- 01 

3.1782 • 

IO- 01 

3.0- 

IO- 01 

1.6872 • 

10-oi 

3.9367 • 

IO- 01 

1.5000 ■ 

IO- 01 

4.0- 

10-oi 

5.3029 - 

Iq-02 

2.5572 • 

IO- 01 

5.9770 ■ 

Iq-02 

5.0- 

IO- 01 

1.3363 • 

Iq-02 

1.6399 • 

IO " 01 

1.9252 • 

l 0 -° 2 

6.0- 

IO” 01 

2.9173 • 

Iq-OS 

9.8636 • 

l 0 -° 2 

4.6135 • 

IQ-OS 

7.0- 

IQ- 01 

6.6102 • 

10 04 

5.2156- 

l 0 -° 2 

6.9805 • 

l 0 -° 4 

8.0- 

io- 01 

1.4527- 

10 04 

2.1510 ■ 

l 0 -° 2 

4.6040 ■ 

10-os 

9.0- 

IQ-Oi 

1.5065 • 

10-os 

4.8033 • 

IQ-OS 

4.1953- 

10-07 


Table 33: Same as table 25: FFNS, ^l 2 f = 10 3 GeV 2 , n 2 R /n 2 F = 2, a s (fi 2 F ) = 0.1188431. 


- 23 - 



































X 

Df(x 

, Go) 

Df(x 

,Go) 

Df{x 

,Go) 

1.0- 

l 0 -°2 

1.0772 • 

lO+oo 

6.0595 • 

10-oi 

-1.8519 

• io- 01 

5.0- 

Iq-02 

1.5140 ■ 

lO+oo 

1.4696 ■ 

lO+oo 

8.7330 • 

10-oi 

1.0- 

10-oi 

9.3978 • 

10-oi 

1.1033 ■ 

lO+oo 

6.6882 • 

10-oi 

2.0- 

to- 01 

3.6193 • 

10-oi 

6.1004 ■ 

10-oi 

3.2441 ■ 

10-oi 

3.0- 

to- 01 

1.2108 • 

io- 01 

3.5792 • 

10-oi 

1.3670 ■ 

10-oi 

4.0- 

to- 01 

3.3810 ■ 

l 0 -° 2 

2.1940 ■ 

10-oi 

4.9203 • 

Iq-02 

5.0- 

to- 01 

8.1014 • 

CO 

o 

1 

O 

1.3386 • 

10-oi 

1.4472 • 

Iq-02 

6.0- 

to- 01 

1.8753 • 

CO 

o 

1 

O 

7.6532 • 

Iq-02 

3.1898- 

l 0 -° 3 

7.0- 

to- 01 

4.6737 • 

i 0 -° 4 

3.8109 ■ 

Iq-02 

4.4482 • 

10” 04 

8.0- 

10-oi 

9.9960 • 

10-° 5 

1.4500 ■ 

Iq-02 

2.7031 ■ 

IQ-OS 

9.0- 

IQ" 01 

8.8903 • 

10-06 

2.8360 • 

l 0 -° 3 

2.5268 • 

10-07 


Table 34: Same as table 25: FFNS, / l 2 f = 10 4 GeV 2 , n 2 R //j' 2 f = 1/2, a s {^? F ) = 0.1096856. 


X 

Dg + (x,Qo) 

Df{x 

3 Go) 

Df{x 

, Go) 

1.0- 

Iq-02 

2.1401 • 

10+°° 

1.5232 • 

10+oo 

7.3458 • 

10-oi 

5.0- 

Iq-02 

1.4826 ■ 

lO+oo 

1.4896 ■ 

10+oo 

8.9453 • 

IO - 01 

1.0- 

IO - 01 

8.9001 ■ 

10-oi 

1.0615 ■ 

10+oo 

6.2736 • 

IO - 01 

2.0- 

10- Ql 

3.4683 • 

10-oi 

5.8798 • 

IO - 01 

3.0245 • 

IO - 01 

3.0- 

IO - 01 

1.1850 ■ 

10-oi 

3.5255 • 

IO - 01 

1.3073 • 

IO - 01 

4.0- 

IQ-Oi 

3.3718- 

Iq-02 

2.2005 • 

10-oi 

4.8411 ■ 

Iq-02 

5.0- 

IO - 01 

8.1802 • 

Iq-OS 

1.3584 • 

IO - 01 

1.4607 - 

Iq-02 

6.0- 

IO” 01 

1.8956 • 

IQ-OS 

7.8336 • 

l 0 -° 2 

3.2876 • 

Iq-OS 

7.0- 

IQ - 01 

4.7094 - 

10 04 

3.9357 • 

l 0 -° 2 

4.6589 • 

l 0 -° 4 

8.0- 

io- 01 

1.0137- 

10 04 

1.5152 • 

l 0 -° 2 

2.8525 • 

10-os 

9.0- 

IQ - 01 

9.1822 • 

Iq-06 

3.0223 • 

CO 

o 

1 

O 

2.4445 • 

Iq-07 


Table 35: Same as table 25: FFNS, /j 2 f = 10 4 GeV 2 , ji 2 R //j? f = 1, a s {^ 2 F ) = 0.1036710. 


X 

Dg ( X,Qo ) 

DZ"(x 

3 Go) 

Df{x 

, Go) 

1.0- 

Iq-02 

2.4103 • 

lO+oo 

1.8049 • 

10+oo 

1.0163 • 

lO+oo 

5.0- 

Iq-02 

1.4223 • 

10+oo 

1.4575 • 

lO+oo 

8.5959 • 

IO- 01 

1.0- 

IQ-Oi 

8.5836 • 

IQ-Oi 

1.0315 • 

lO+oo 

5.9413- 

10-oi 

2.0- 

10- Ql 

3.4550 ■ 

10-oi 

5.8112 • 

IO- 01 

2.9304 - 

IO- 01 

3.0- 

IO- 01 

1.2159 - 

10-oi 

3.5500 - 

IO- 01 

1.3022 • 

IO- 01 

4.0- 

10- Ql 

3.5427 • 

Iq-02 

2.2467 • 

IO- 01 

4.9424 • 

Iq-02 

5.0- 

IO- 01 

8.7151 ■ 

IQ-OS 

1.4012 • 

IO- 01 

1.5228 • 

Iq-02 

6.0- 

IO” 01 

2.0155 • 

IQ-OS 

8.1569- 

l 0 -° 2 

3.4896 • 

l 0 -° 3 

7.0- 

IQ- 01 

4.9631 - 

Iq-04 

4.1424 • 

l 0 -° 2 

5.0279 - 

l 0 -° 4 

8.0- 

io- 01 

1.0723 • 

10 04 

1.6180- 

l 0 -° 2 

3.1295 • 

10-os 

9.0- 

IQ- 01 

9.8900 ■ 

10 06 

3.3055 • 

IQ-OS 

2.6762 - 

10-07 


Table 36: Same as table 25: FFNS, ^l 2 f = 10 4 GeV 2 , [i 2 R /fi 2 F = 2, a s (fi 2 F ) = 0.0982996. 
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This script will automatically produce the relevant numbers and plots. Note that a 
previous installation of the APFEL library version 2.3.0 or higher and the Gnuplot plotting 
tool is required. For the installation of APFEL please refer to: 

http://apfel.hepforge.org/download.html 
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